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ABSTRACT 

a^: 

^ I We present a comparison of semi-analytic models of the phase-space structure 

_H ! of tidal debris with measurements of average distances, velocities and surface 

^ I densities of stars associated with the Sagittarius dwarf galaxy (Sgr), compiled 

I from all observations reported since its discovery in 1994. 
CN I We find that several interesting features in the data can be explained by 

I I these models. The properties of stars ~ ±10 — 15° away from the center of Sgr 

K* ' — in particular, the orientation of material perpendicular to Sgr's orbit (Alard 

if^ [ 1996, Mateo et al 1996, Fahlman et al 1996, Alcock et al 1997) and the kink in 

^ I the velocity gradient (Ibata et al 1997) — are consistent with those expected for 

O ! unbound material stripped during the most recent pericentric passage ~ 50Myrs 

2? I ago. The break in the slope of the surface density seen by Mateo, Olszewski 

"^ I & Morrison (1998) at 6 ~ —35° can be understood as marking the end of this 

Qh| material. However, the detections beyond this point are unlikely to represent 

O ' debris in a trailing streamer, torn from Sgr during the immediately preceding 

'^ I passage ~ 0.7Gyrs ago, as the surface-density of this streamer would be too 

f^ I low compared to observations in these regions. The low b detections are more 

. ^ ' plausibly explained by a leading streamer of material that was lost more that 

/\ ' 1 Gyr ago and has wrapped all the way around the Galaxy to intercept the 

c^ ■ line of sight. The distance and velocity measurements at 6 = —40° reported in 
Majewski et al (1999) also support this hypothesis. 
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We determine debris models with these properties on orbits that are 
consistent with the currently known positions and velocities of Sgr in Galactic 
potentials with halo components that have circular velocities fcirc = 140 — 200 
km/s. In all cases, the orbits oscillate between ~12 kpc and ~40 kpc from the 
Galactic center with radial time periods of 0.55 — 0.75Gyrs. The best match 
to the data is obtained in models where Sgr currently has a mass of ~ IO^Mq 
and has orbited the Galaxy for at least the last 1 Gyr, during which time it has 
reduced its mass by a factor of 2-3, or luminosity by an amount equivalent to 
~ 10% of the total luminosity of the Galactic halo. These numbers suggest that 
Sgr is rapidly disrupting and unlikely to survive beyond a few more pericentric 
passages. These conclusions are only tentative because they rely heavily on 
the less-certain measurements of debris properties far from the center of Sgr. 
However, they demonstrate the immense potential for using debris to determine 
Sgr's dynamical history in great detail. 



Subject headings: galaxies: individual (Sagittarius) - galaxies: dwarf - galaxies: 
interactions - galaxies: kinematics and dynamics - Galaxy: kinematics and 
dynamics - Galaxy: evolution - methods: analytical 



Introduction 



Since its discovery in 1994 (Ibata, Gilmore & Irwin 1994, 1995 — hereafter [I94| and 



1951) , the history of the Sagittarius dwarf galaxy (Sgr) has been a matter of some debate. 
Sgr is centered on the globular cluster M54 at (/, b) = (6°, —14°) on the opposite side of 
the bulge from us at 16 kpc from the Galactic center and is apparently moving towards 
the Galactic plane ( |I94| , |95|, Ibata et al 1997 — hereafter ^9^ ). The isopleths revealed 
by careful star count analyses (p.94|, [19 5|, |197|) are highly elongated along the direction of 



its motion, and populations plausibly associated with it have been found even further 
upstream and downstream (see Alard 1996, Mateo et al 1996, Fahlman et al 1996, Alcock 
et al 1997, Mateo, Olszewski & Morrison 1998 and Majewski et al 1999 — hereafter 
A96t ^9|, H, ^, ^9| and Paper H respectively). In addition, [Irwin fc Totten (1999) 



have found an excess in the surface density of carbon stars aligned with the direction of 
elongation of Sgr that completely encircles the Galaxy. This suggests that Sgr is being 
strongly influenced by the Milky Way's tidal field, possibly completely destroyed, yet its 
proximity to the Galactic center also indicates that it may have had time to circle the 
Galaxy many times. If we are currently witnessing just one in a long series of encounters it 
seems unlikely that it should also be the last. 
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Johnston, Spergel & Hernquist (1995) tackled this problem using A^-body simulations 
to investigate the evolution of stellar systems along highly eccentric orbits that would suffer 
few encounters during the lifetime of the Galaxy. They found that they could fit nearly all 
the observations available at that time with a wide variety of models, without requiring 
that this be either the first or last passage of Sgr. However, they did not attempt to fit 
the lack of velocity gradient reported in the center of Sgr ([195|) . [Velasquez fc White (1995 



demonstrated that this datum alone suggests that Sgr cannot be on a highly eccentric orbit, 
and their estimate for the proper motion of Sgr has since been confirmed by observations 
( [19 71) , and rules out the specific models investigated by [Johnston, Spergel fc Hernquist] 



1995)[ . [19 7[ argue that if Sgr is embedded in a dark matter halo, it could in fact be 
impervious to tidal destruction, but concede that it could be being stripped and limited 
by the Galactic tidal field. This idea is illustrated with numerical simulations by [Ibata fc 



[Lewis (1998)[ Finally, [Zhao (1998)[ suggests that the proper motion measurements of Sgr 
and the Large Magellanic Cloud (LMC) are consistent with a scenario in which Sgr was 
initially on an orbit much farther from the Galactic center - and hence immune to the affect 
of the Galactic tidal field - and has only recently been flung into it's current orbit through 
a collision with the LMC. 

If we could measure how fast Sgr is currently losing mass, and constrain how much 
mass it has lost in the past, we could make some progress towards resolving this debate 
about its history and current state. In this paper we address these questions by building 
semi-analytic models of the structure of tidal debris that are consistent with all the current 
position- and velocity-space observations of material believed to be associated with Sgr. In 
§2.1 we describe sources and characteristics of the data used to constrain the models and in 
§2.2 we outline the method for building the semi-analytic models. We present models which 
agree with all the observations in §3 and discuss the extent of the limitations that the data 
places on Sgr's orbit, current mass and mass loss history, and the potential of the Milky 
Way. We summarize our results in §4. 



2. Tools 

2.1. Summary of Data Sets 

In this study we have combined data from eight different investigations that are 
currently available in the literature. We take the distance (24 kpc), line-of-sight velocity 



(170 km/s) and proper motion from [194[ , [19 5| and [19 7[ , and the velocity trends along the 
major axis of Sgr from Table 3 of [19 7[ . At b = —7° and b = —4° (+7° and +10° away from 
the center of Sgr towards the Galactic disk) we have two more distance estimates from the 
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discovery of associated RR Lyraes ( A96 , A97 ), and at 6 = —23° we have distances both from 



RR Lyraes ( |M96| ) and main- sequence fitting (|F96|) to Sgr populations found serendipitously 



in studies of the globular cluster M55. We also have one distance and a tentative velocity 
measurement at 6 = —40° (see [Paper i| for discussion) . Finally, we use the surface density 



profile reported in |M98| , which was produced from star counts at the position of the Sgr 
upper-main sequence in a set of 24 color-magnitude diagrams that cover the expected locus 
of Sgr's debris trail sky from b = —24° to b = —48°. 

The symbols in Figure |l| summarize this combined data set, with the solid line in Panels 
(a), (b) and (c) showing the locus of an orbit calculated from the position and velocity of 
Sgr and the meaning of the symbols in each panel summarized in Panel (c). Panel (a) shows 
the position of Sgr projected onto the plane containing the Sun (at (X, Z) = (—8, 0)kpc 
in this coordinate system) and perpendicular to the Galactic disk. Panels (b), (c) and 
(d) show the distance from the Sun, line-of-sight velocity in a Galactic frame of rest and 
(arbitrarily normalized) surface density measurements plotted against Galactic latitude. 
(Note: the surface density of stars identified as possibly belonging to Sgr in [Paper I| agrees 



with |M98 's estimate and this point is also included in Panel (d) for comparison.) 



When we allowed Sgr's proper motion to vary within the quoted error bars ( [197] ) and 
compared the different data sets both with each other and with the loci of the orbits we 
found a number of interesting features, already noted by previous authors: the distances in 
the central regions of Sgr follow a line perpendicular to the orbit rather than aligned with 
it; there is a kink in the velocities measured by [19 7[ that is hard to reconcile with the orbit 
calculations alone; although the distance measurement at 6 = —40° from Paper I typically 
lies close to the orbit, the velocity measurement does not; and there is a break in the slope 
of the Sgr star counts at 6 = —35°. In §3 we will demonstrate how all these features could 
be a natural consequence of debris dispersal from the disruption of Sgr. 



2.2. Summary of Method 

In this section we briefly outline the semi-analytic method we employ to model debris 
dispersal — the details are reported in [Johnston (1998a, 1998b)[ The method predicts 



the full phase-space positions and densities of mass lost at a given time from a satellite of 
known mass, rrisat, and motion orbiting in a specified parent galaxy potential, $, at any 
later time. 

When a satellite disrupts, escaping stars that lose (gain) energy move ahead of (behind) 
the satellite to form leading (trailing) streams of debris along its orbit. In a given parent 
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galaxy potential $, the gain or loss of energy AE" of a star can be related to its position 
at any time following the mass loss event by comparing the time period of a circular orbit 
of energy E + AE with that of the unperturbed satellite orbit of energy E (ignoring the 
weak dependence of orbital periods on angular momentum), and adjusting the azimuthal 
position of debris along the orbit accordingly. The star's AE also gives an estimate of its 
offset in distance from the orbit: AR ~ AE / [d^ / dR) . The semi-analytic method combines 
this approximate modeling of particle dynamics with (1) an analytic description of the 
energy scale e = R{msa.t/ McaiY^^d^ / dR (where Moai is the mass of the parent galaxy 
enclosed within radius R) in tidal debris; (2) the intrinsic energy distribution (in scaled 
energy units AE/t) observed in A^-body simulations of satellite disruption. The method is 
further simplified by the observation that in simulations of satellite destruction most mass 
loss occurs near the pericenter of the satellite's orbit and hence its mass loss history can 
be approximated as a series of discrete mass loss events. Despite these simplifications, the 
semi-analytic method successfully reproduces the density, position and velocity of streamers 
in a wide variety of simulations (for a given satellite mass, orbit and approximate mass-loss 
history) at a tiny fraction of the simulation's computational cost (Johnston 1998a, 1998b). 
This makes the method an ideal tool for interpreting observations such as those described 
above because it allows a quick and thorough exploration of parameter space. 

Several limitations should be noted. First, the method does not take into account the 
effect of dynamical friction on the orbit of the satellite. From [Binney fc 'I'remaine (198T)] 



equation 7-27, the dynamical friction timescale for a IO^Mq satellite orbiting at i? =15 kpc 
from the center of the Galaxy is ~ 5 Gyrs. Since we consider < 1.5 Gyrs in the evolution of 
satellites of mass < lO^M© on orbits where most time is spent at i? > 15 kpc, we expect this 
simplification to be unimportant. Moreover, our interpretation of the data is only intended 
to be approximate. However, once our knowledge of the structure and extent of debris 
associated with Sgr becomes detailed enough to warrant a more specific interpretation of its 
history, the effect of dynamical friction should be included. 

Second, the method is most accurate for small satellites, with masses such that 
{'msa.t/MGaiY^'^ <^ 1 — for our own Galaxy this suggests it should only be used to study 
satellites, such as Sgr, with msat < lO^M©. Finally, the method was developed and tested 
under circumstances where the satellite's orbit was eccentric but not radial, and where the 
satellite's structure was supported by random motions. Sgr satisfies both these conditions. 



Results 
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3.1. General Approach 

Our aim is to explore whether the combined data sets are consistent with the 
semi-analytic models of debris dispersal from satellite disruption and how much such a 
model can tell us about Sgr's history, orbit and the Milky Way's potential. We do not 
attempt to use an automated routine to find the "best fit" model since we believe that 
neither the extent of the data nor the accuracy of the models at present justify this level of 
sophistication. We instead simply compare the data and the models by visual inspection, 
varying four parameters: 

1. / — the fraction of mass lost on each pericentric passage. We assume that mass loss 
occurs instantaneously at pericenter and that / is constant (as is approximately true 
in simulations). 

2. TTisgr — the mass of material still bound to Sgr and centered around M54. 

3. ft — Sgr's velocity transverse to the line of sight. We take Sgr's orbit to be exactly 
polar, fix its line-of-sight velocity to its observed value (170 km/s), and allow Vt to 
vary within the error-bars on the proper motion measurement quoted in 197. In 
fact, Sgr's motion parallel to the Galactic plane is poorly constrained by current 
observations ( [19 7| ), although the surface density contours do indicate that it must be 
on a near-polar orbit, and hence it is likely to be small. We will discuss the effect of 
non-polar orbits in flattened halo potentials in §3.5. 

4. fcirc — the circular velocity of the dark matter halo, as defined in equation (|l|) below. 
We integrate Sgr's orbit in Milky Way potentials of the form $ = $buige + "^disk + ^haio 
where 



•^ bulge 
$disk 



GM, 



bulge 



VWT^ + c 



R^ + {a + V^^+¥) 



2 



,2 



$halo = ^HR' + z^q' + d'). (1) 

Here R and z are cylindrical coordinates aligned with the Galactic disk and the 
numerical values adopted for the parameters are Mdisk = 1-0 x 10^^,Mbuige = 
3.4 X 10^°, a = 6.5, b = 0.26, c = 0.7 and d = 12.0, where masses are in Mq, velocities 
are in km/s and lengths are in kpc. We vary fcirc between 140 km/s and 200 km/s 
and set g = 1 to consider spherical halo potentials. The effect of flattened halo 
potentials is discussed in §3.5. The top panel of Figure ^ plots the rotation curves for 
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the combined potential $ with fcirc = 140, 170 and 200 km/s, and the bottom panel 
shows the mass enclosed for these choices within a given radius. 

In Figure ^ we illustrate the general characteristics of debris dispersal by 
following a semi-analytic model for the evolution of a satellite with parameters 
(/, msgr, ft, fcirc) = (0.3, IO^Mq, 242km/s, 174km/s) over three pericentric passages. Each 
panel is centered on the Galaxy and 100 kpc on a side. The solid line shows the past 
and future path of the satellite, and the hexagon shows its position. The left hand panels 
correspond to instants shortly after each pericentric passage (with the bottom left hand 
panel corresponding to Sgr's current position and velocity) and the panels in each row are 
equally spaced between these points in time. The grey ribbons show the location of the 
debris lost on each pericentric passage. The darkest and shortest pair of ribbons in the 
final panel correspond to the leading and trailing streamers from the most recent passage 
(hereafter referred to as passage or by ripass = 0) and the paired streamers from the 
preceding passages ( hereafter referred to by ripass = —1,-2...) are progressively lighter and 
longer. Note that the large msgr means that the tidal scale (msgr/Mcai)^''^ ~ 1/10 is large 
and so the leading (trailing) streamers are on orbits that come (go) considerably closer to 
(further from) the Galactic center than the satellite itself. Each of these streamers will have 
a width less than, but comparable to the amount by which it is offset from the orbit so there 
will considerable overlap between the populations whose average positions are indicated by 
the grey lines. 

In Figure ^ we repeat Figure |^, but with the predictions of the semi-analytic model 
shown in Figure ^ overlaid. In the following subsections we use this to illustrate the general 
successes of using the dynamics of tidal debris to explain the data sets (§3.2). We then go 
on to discuss whether the current observations can place any stringent limits on the varied 
parameters (§3.3), and outline the implications for Sgr's mass loss history (§3.4). Finally, 
we examine to what extent our discussion will remain valid for non-polar orbits in flattened 
halo potentials (§3.5). 



3.2. Successes of the Models 

The location of the streamers (grey ribbons) and the orbit (solid line) projected onto 
the plane containing the Sun and perpendicular to the Galactic disk are shown in panel (a) 
of Figure |[ Note that the streamers from the most recent pericentric passage (npass = 0) 
are nearly perpendicular to the orbital path, and that those from the preceding passages 
(^pass = — 1, —2...) are clearly separated by a noticeable gap ahead and behind the satellite 
along its path. This general configuration is apparent in all models, and results from Sgr 



being only slightly past pericenter in orbital phase — once it reaches apocenter the streamer 
from passage will spread far enough to overlap in phase with the older ones and the debris 
will be continuous (see Figure ^. However, currently this means that as we trace the debris 
back along the orbit of Sgr beyond the extent of material unbound during passage we 
will not immediately detect material from passage —1, but rather debris that has wrapped 
entirely around the Galaxy and hence must have been torn from Sgr on passages with 
^pass < —2, as illustrated by the position of the lightest leading streamer in the Figure. 

In panel (b) of Figure ^ we plot the distances to the streamers as seen from the Sun in 
the direction of the center of the Galaxy as a function of Galactic latitude. Note that the 
streamers from passage agree very well with the peculiar orientation of material reported 
in [A96| , [[V196| , |1^'96| and |A97| . The distance estimate we found for our field further down the 



streamer in Paper I agrees with the prediction for material in the leading streamer from 
passage —2 that has wrapped entirely around the Galaxy to intersect our line of sight. 
This point serves as a discriminant between the streamers: trailing streamers are offset 
from the satellite's path to greater distances and should not be observed so close to the 
Galactic center; and material in leading streamers from passages with npass ^ ~3 will lie 
progressively closer to the satellite's orbit (and further from us along the line of sight close 
to Sgr) at a given orbital phase. These arguments apply to all models. 

In panel (c) of Figure ^ we plot the line-of-sight velocity as seen from the Sun in the 
Galactic rest-frame. The model agrees with the kink in the velocity gradient for the points 
close to the center of Sgr that was commented on by |I97| . The velocity reported at b = —40° 
in Paper I is too low to agree with a direct extrapolation of the orbit of Sgr over 20° within 
the stated error-bars of the proper motion in any of our Galactic models. However, it can 
be explained consistently within the model where this material is part of a leading streamer 
from passage —2. For a given Galactic model, the distance and velocity of this point, 
provides the strongest constraint on the orbit of Sgr. 

In panel (d) of Figure § we show the surface density of each streamer, normalized so 
that the streamer from passage has the same amplitude as the data in the region reported 
in M98. In all models, there is a gap of 10° — 15° between the leading and trailing streamers 
from passage corresponding to the location of mass still bound to Sgr. The models suggest 
that the break in slope of the surface density at 6 ^ —35° corresponds to the end of the 
streamer from passage 0. The location of this break is dependent on rrisgr and the density 
of debris beyond this point is sensitive to /. 
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3.3. Limits on Parameters 

In order to investigate to wliat extent tlie data could limit the parameters we looked 
at a set of potentials with halo components having 140km/s < f drc < 200 km/s. In each 
potential we found the value of Vt that allowed the closest fit (by visual inspection) to all 
the distance and velocity data, except for the tentative velocity measurement at 6 = —40° 
(see 



Paper I| for discussion of this point). We considered debris from the last 3 passages 
alone and did not look at debris from earlier passages because (as noted in the previous 
section) these streamers would be too distant to fit the point measured at 6 = —40° in 
[Paper 1| . In all potentials we could find a value of Vt that provided a reasonable fit to these 
data within the error bars of the proper motion measurement reported in 197. As f drc was 
increased from 140 km/s to 200 km/s, the necessary vt increased from 205-275 km/s (with a 
range of about ±10 km/s about each value) and the time periods and time interval since the 
most recent pericentric passage decreased from 0.75-0.55 Gyrs and 60-40 Myrs respectively. 
However, the pericenters and apocenters of the orbits orbits lay in narrow ranges (12-14 kpc 
and 40-43 kpc). We conclude that although the current observations do not significantly 
refine our model of the Galactic potential, the proper motion measurement of Sgr, nor the 
prediction for the time period of its orbit, they do constrain the orbital shape. If we knew 
the proper motion of Sgr with greater precision we could use the properties of the debris to 
constrain the Galactic potential more tightly. 

Note that the model that fits all the distance and [I97f s velocity measurements 



also provides a good match to the uncertain velocity measurement reported in Paper I| , 
even though it was not considered as a constraint. This success supports the tentative 
identification in [Paper I| of a group of stars moving at this velocity as indeed being 
associated with Sgr. 

Once we determine a value of Vt in each potential we can use the surface density 
measured in M98 to constrain rrisgr- Increasing msgr causes the length of each streamer to 
increase (approximately oc rnj^j.). In all potentials considered we found that msgr ~ lO^M© 
in order for the streamer from passage to extend to the point where the slope in the 
surface density profile in M98 changed. M98 estimate a total luminosity associated with Sgr 
of < 8 X IO^Lq (i.e. including their detections) so this suggests an overall mass-to-light ratio 
of > 10, which is consistent with observations of other dwarf spheroidal satellites in orbit 
around the Milky Way, and with the independent estimates for its mass by 197. However, 
this large mass is somewhat inconsistent with the measured velocity dispersion of ~ 10 
km/s along the entire length of Sgr ( [197[ ). This inconsistency could arise because the finite 
size of the satellite is ignored in the semi-analytic models, where the debris is effectively 
being released from the satellite's center. This simplification would affect the predictions 
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for the length of the tidal streamers from the most recent pericentric passage (passage 0) 
whose size is comparable to that of the satellite itself, but should be unimportant for the 
much longer streamers from earlier passages. 

Finally, we found that / was required to be in the range 0.2 < / < 0.3 in order for the 
relative amplitude of debris at either side of the break in the slope of the surface density 
profile approximately to be fit by material in the streamers from passages and —2. 



3.4. Mass Loss History and Future Survival 

One question that has been addressed repeatedly since Sgr's discovery is whether we 
are witnessing its death throes or will Sgr survive another passage. The models that best 
fit the data described in the previous section had 0.2 < / < 0.3, which means that Sgr 
has lost between half and two-thirds of its mass over the last three pericentric passages (or 
1.1-1.5 Gyrs depending on the Galactic potential). We can crudely check this picture by 
assuming that our success in modeling the distance and velocity features means that Sgr is 
dominated by unbound material beyond at least 10° from its center. M98 report that the 
surface brightness (in mag/arcsec~^) of their points between 10° and 20° from Sgr can be 
fitted by a function S^ oc 0.23^ (where 9 is the separation from the center of Sgr in degrees) 
which agrees with the central surface brightness of Sgr when extrapolated inwards. If 
mass follows light, this suggests that the mass surface density will fall as S oc exp(— ^/lO). 
Hence, if Sgr and its debris have roughly constant width along its orbit the ratio of the 
mass between 10° < 9 < 20° (i.e. in the streamers from passage 0) to that within 9 = 10° 
will be ~ exp{—l) = 0.37, in rough agreement with our assumed fractional mass loss rates. 
Such a large mass loss rate rules out the picture where Sgr is merely distorted by the Milky 
Way's tidal field with very little evolution and suggests that Sgr is unlikely to survive many 
more pericentric passages. This mass loss rate corresponds to a decrease in Sgr's luminosity 
in the last 1.1-1.5 Gyrs alone that is of order 10% of the total luminosity of the stellar halo 
(~ IO^Lq, see freeman 1996|) . 



3.5. Halo Flattening and Non-Polar Orbits 

We have so far limited our discussion to polar orbits in potentials with exactly spherical 
halos. Non-polar orbits in a flattened halo potential (i.e. setting g < 1 in equation |l|]) will 
precess and disturb the alignment of the leading and trailing streamers. In principle, we 
could use the apparent precession of the debris at different points along the streamers to 
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place limits on both the direction of the proper motion and the degree of flattening of the 
Galaxy's potential. However, the strongest constraints on these quantities would come from 
detections (or lack of detections) farthest from the center of Sgr (e.g. |M98| , [Paper ij ) which 



are currently the most uncertain. Hence, we do not try to place limits on these dimensions 
of parameter-space, and instead simply examine to what extent our discussion of debris 
evolution in a spherical halo will remain valid in the non-spherical case. 

Figure | contours the projected angular difference between Sgr and a point 27r ahead 
along it's orbit. The orbits were run in halo potentials with different q and with initial 
radial velocity and motion perpendicular to the Galactic disk set by [197f s observations, but 
with motions perpendicular to the line of sight and parallel to the Galactic disk v^ allowed 
to vary. The positive (solid) and negative (dotted) contours are spaced by 5 degrees. This 
plot provides upper limits to the apparent orbital precession of material in the leading 
streamer since debris on orbits more tightly bound to the Galaxy than Sgr would lie closer 
to the Galactic center than it's orbit and hence have smaller angular separations from Sgr 
when seen from the Sun in projection. Figure § clearly cautions that our identification of 
M9S 's and Paper l| 's detections as material in the leading streamer could be wrong if g < 0.9 



and If^l > 20 km/s. Although this condition appears to be rather restrictive, note that this 
modest potential flattening translates to an axis ratio ~ 0.7 in the corresponding density 
distribution (see Binney & Tremaine 1987, equation 2-55b). In comparison, a recent study 
combining limits from stellar-kinematics and from the flaring of gas out of the Galactic 
plane constrains the dark matter halo to have an axis ratio in density of 0.75 ± 0.25 ( piling 
|fc Merrifield 1998|) . 



4. Summary and Conclusions 

In this paper, we have combined all observations of Sgr that have been reported in the 
literature since 1994 and directly compared them to semi-analytic models of the dispersal of 
material once it becomes unbound from a satellite. By varying the circular velocity of the 
halo component of our Milky Way potential along with Sgr's transverse velocity, current 
mass and mass loss rate, we have found that: 

1. The current data cannot constrain the Milky Way potential and Sgr's transverse 
velocity independently, but do favor specific orbital characteristics: the orbits of the 
best fit models all had pericenters ~ 13kpc, apocenters ~ 41kpc, radial time periods 
~ 0.55 — 0.75Gyrs and time since the most recent pericenter ~ 40 — 60 Myrs. If the 
proper motion of Sgr were known with greater precision, the data could be used to 
constrain the Galactic potential more tightly. Note that [Murali fc Dubinski (1999 
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also found that poor knowledge of proper motion was the hmiting factor to using 
globular cluster tidal tails to measure the Galactic potential. 

2. All the peculiar features in the data could be explained with the dynamics of unbound 
material: 

• The orientation of material perpendicular to Sgr's orbit close to its center 
reported in [A96| , |A97| and Paper 1| , and the kink in the radial velocities reported 



in P97| are consistent with the predicted locus of debris from the most recent 
pericentric passage. 

The radial velocity and distance measurements of material at 6 = —40° reported 
in Paper I correspond to those expected for material in a leading streamer 
that has wrapped entirely around the Galaxy to intercept our line of sight. 
The density of this material in this region dominates over trailing material 
because there is a clear gap between material stripped from the most recent and 
immediately preceding pericentric passages. 

The break in the slope of the surface density profile reported in [M98| naturally 



occurs in our model at the transition between material in the trailing streamer 
from the most recent passage and the leading streamer that has wrapped around 
the Galaxy. 

3. Models that are consistent with the data have current masses of Sgr of ~ IO^Mq, 
and mass loss rates such that Sgr was 2-3 times larger < l.SGyrs ago. This mass 
loss rate corresponds to a decrease in Sgr's luminosity that is of order 10% of the 
total luminosity of the stellar halo and suggests Sgr is unlikely to survive many more 
pericentric passages. 

We consider the general scenario outlined in our first two points to be fairly robust, 
provided the halo potential is not highly flattened (see §3.5). The specific model for the 
mass and mass loss rate in the third point is more tentative. Future observations should 
reveal the history of Sgr in greater detail. Velocity measurements along the entire length 
of [lV198|'s data set would clearly confirm or rule out our identification of this material as 



part of the leading streamer and distance estimates would indicate whether we are indeed 
seeing only material wrapped around the Galaxy from the pericentric passage ~ 1.1 — 1.5 
Gyrs ago (i.e. npass = —2), or if we might also be seeing material (a few kpc more distant) 
from earlier passages. The absence of debris from these earlier passages in M98's data set 
could provide support for Zhao's (1998) model of Sgr's history, which requires Sgr to be on 
a much less tightly bound orbit before suffering a collision with the LMC 2 Gyrs ago. This 
scenario limits the number of close encounters with the Milky Way that Sgr has undergone. 
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and hence the number of streamers from different passages associated with it. On the 
other hand, if Zhao's model can clearly be ruled out, we instead face the task of building a 
consistent picture in which Sgr has orbited the Galaxy many times and is only now being 
completely destroyed. 

Additional detections (ideally with distance, velocity and surface density estimates) or 
lack of detections of stars associated with Sgr in fields in the Northern Galactic hemisphere 
(e.g. [Paper 1| ) or radial velocity measurements along |M98| 's data set would allow us to 



pin down Sgr's path, which in turn would more precisely constrain both the shape and 
depth of the Galactic potential and Sgr's proper motion. Unfortunately, the converse of 
this statement is that our poor knowledge of Sgr's motion parallel to the Galactic disk 
and our uncertainty about the degree to which the Galaxy's potential is flattened mean 
that any predictions of where to find such streamers are very imprecise. However, the 
recently discovered carbon star stream (Irwin & Totten 1999) could provide a starting point 
for such studies — a preliminary analysis of the carbon star data set alone suggests that 
the halo potential is near spherical (private communication from Rodrigo Ibata). Further 
into the future the GAIA satellite (Global Astrometric Imager for Astrophysics — to be 
launched in 2009) will provide proper motions of all stars brighter than 20th magnitude 
with ~ 5 — 200/ias/yr accuracy and allow us easily to distinguish Sgr's debris across the 
sky from foreground contaminants. 

Finally, we caution that the results reported here are only intended as a guide to the 
probable history of Sgr and that more detailed velocity and distance information along 
its streamers is needed to come to firmer conclusions. The semi-analytic models are not 
exact, and, with more data to use as constraints, any results should ideally be confirmed 
using N-body simulations that also take into account the effect of dynamical friction. 
Nevertheless, this study illustrates the tremendous power of using debris, rather than the 
properties of bound material, to understand the history of satellite galaxies. 
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from the Institute for Advanced Study. SRM is supported by a Cottrell Scholar Award 
from The Research Corporation, and a David and Lucile Packard Foundation Fellowship. 
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Fig. 1. — Summary of data used in this investigation. Panel (a) shows the current position 
(triangle) and a possible orbit (sohd line) of Sgr projected onto the plane perpendicular to 
the Galactic disk and containing the Sun at (X, Z) = (—8, 0)kpc. The arrow indicates its 
direction of motion. Panels (b), (c) and (d) show the distance from the Sun, line-of-sight 
velocity in a Galactic rest frame and surface density of the data points plotted against their 
Galactic latitude. The orbit is again overlaid in Panels (b) and (c). The sources of the data 
points are coded with the same symbols in each panel, as labeled in Panel (c). 
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Fig. 2. — Top panel — the velocity of a circular orbit in the disk plane at distance R from 
the Galactic center for the Galactic model given in equation (Q) with fcirc = 140, 174, 200 
km/s. Bottom panel — the mass enclosed within distance R for the same models. 
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Fig. 3. — Evolution of debris lost from a satellite over three pericentric passages projected 
onto the plane perpendicular to the Galactic disk and containing the Sun. In each panel the 
solid line shows the past and future path of the satellite, and the hexagon shows its position. 
The panels are equally spaced by 0.22 Gyrs in time, centered on the Galaxy and lOOkpc 
on a side. The grey lines show the positions of the streamers. The bottom left hand panel 
corresponds to the position and velocity of Sgr today. 
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Fig. 4. — As Figure |l|, but with the semi-analytic model also overlaid. Each pair of grey 
lines represent the locus of a pair of tidal streamers, with the darkest pair being torn from 
the satellite on the most recent pericentric passage, 50 Myrs ago. 
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Fig. 5. — Positive (solid) and negative (dotted) contours of angular precession of orbit, as 
a function of the flattening q of the halo potential and the velocity of Sgr perpendicular to 
the line of sight and parallel to the plane of the Galactic disk f^. The levels are spaced by 
5 degrees. 



